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In this work, we investigate the classical loop models doped with monomers and dimers on a square lattice, 
whose partition function can be expressed as a tensor network (TN). In the thermodynamic limit, we use the 
boundary matrix product state technique to contract the partition function TN, and determine the thermodynamic 
properties with high accuracy. In this monomer-dimer-loop model, we find a second-order phase transition 
between a trivial monomer-condensation and a loop-condensation (LC) phases, which can not be distinguished 
by any local order parameter, while nevertheless the two phases have distinct topological properties. In the 
LC phase, we find two degenerate dominating eigenvalues in the transfer-matrix spectrum, as well as a non¬ 
vanishing (nonlocal) string order parameter, both of which identify the topological ergodicity breaking in the 
LC phase and can serve as the order parameter for detecting the phase transitions. 

PACS numbers: 64.60.Cn, 05.50.+q, 05.10.Cc, 64.60.F- 


Introduction .— Two dimensional (2D) monomer-dimer 
model has a quite venerable history in statistical mechanics 
[1—4]. The monomer-dimer model can be used to describe the 
absorption of molecules on the surface: the molecule can oc¬ 
cupy two nearest neighboring sites and form a dimer, while 
the empty site is regarded as a monomer [1]. The monomer- 
dimer model can also be related to other statistical models like 
Ising and height models [2, 5], etc, thus it plays the role as a 
quite fundamental statistical model. On a square lattice, the 
fully packed dimer model is found to possess algebraic decay¬ 
ing dimer-dimer correlation, however, doping the system with 
monomers will drive the system out of the criticality and no 
phase transition occurs in a non-interacting monomer-dimer 
model at finite temperatures [6-8]. On the other hand, if one 
introduces pairing interactions between the dimers, there exist 
phase transitions between the low-7’ ordered phase and high- 
T disordered one (Kosterlitz-Thouless type for fully packed 
case [6], and second-order after monomer doping [7, 8]). 

Loop models are also widely studied in statistical mechan¬ 
ics, which is relevant for realistic physical systems and also 
constitutes a quite fundamental mathematical problem [9, 10]. 
The loop structure also plays an important role in certain 
quantum cases, like in the ground state of toric code [11], the 
string-net model [12], and the resonating Affleck-Kennedy- 
Lieb-Tasaki loop spin liquid states [13], etc. In Ref. 14, 
Castelnovo and Chamon couple the toric code model to a ther¬ 
mal bath, and consider the thermal superposition of all possi¬ 
ble loop coverings. They found that the concept of topologi¬ 
cal order also applies in this classical loop system, where the 
low-energy phase space decompose into several distinct topo¬ 
logical sectors. The existence of distinct topological sectors 
breaks ergodicity. One needs to create/annihilate a loop with 
length propositional to system size, which has huge energy 
cost and rare probability to happen, in order to tunnel from 
one sector to another, it thus leads to the topological glass 
behavior [15]. The notion of topological entropy can be gen¬ 
eralized to detect such nontrivial topological order in classi¬ 


cal systems, by noticing that the topological constraint would 
also reduce the entropy in the classical case [14]. Recently, 
Hermanns and Trebst have generalized this entropy character¬ 
ization to general classical string-nets and verify that there are 
corresponding universal topological corrections in the Renyi 
entropy for a number of SU(N)* anyonic theories [16]. 

In this work, we combine the two classical models and in¬ 
troduce a monomer-dimer-loop (MDL) model on a square lat¬ 
tice. The MDL model has a rather compact tensor network 
(TN) representation with a small bond dimension (D = 3), 
and is thus amenable to high precision TN numerical simula¬ 
tions. TN-based numerical methods have been widely used to 
tackle statistical models and have been proved to be a very ac¬ 
curate and reliable tool [8, 17-19]. Through the TN numerical 
simulations, we show that the MDL model has a trivial dis¬ 
ordered phase and a topologically ordered loop-condensation 
(LC) phase, with a second-order transition separating them. 
In addition, we characterize the LC phase with the vanishing 
gap of transfer-matrix spectrum and a nonlocal string order 
parameter (SOP), both of which can be used to pinpoint the 
phase transition. 

Model and method .— Snapshots of several classical con¬ 
figurations in different phases of MDL model are shown in 
Fig. 1. Summing over all possible classical configurations, 
we have the partition function 

H = X exp [~f3(pN m + vN b + uN d )], (1) 

M 

where [c] means the set of all classical monomer-dimer-loop 
configurations, N m is the total number of vertices occupied by 
a monomer, which has an energy of /a; Nb counts the number 
of edges occupied by a loop, and v is the energy per bond of 
a loop; N d is the total number of vertices linked by a dimer 
(with energy per dimer as 2m). In the following, v = 1 is set as 
the energy scale if not otherwise specified. 

The partition function of the monomer-loop model has a 
simple TN representation which form a 7r/4 tilted square 
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FIG. 1: (Color online) Snapshots of the classical phases of the MDL 
model on a square lattice: (a) The trivial monomer-condensation, 
(red) dots are the monomers, A.B label the two sub-lattices; (b, 
c) Loop condensation in the monomer-loop and monomer-dimer- 
loop cases, respectively; (d) Loop condensation in the branching 
monomer-loop model (i.e., classical string-nets). 
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FIG. 2: (Color online) (a) TN representation of the partition function, 
on a tt/ 4 tilted square lattice, (b) The construction of the transfer 
matrix M on a cylindrical geometry by contracting a column of (w) 
rank-4 tensors. The dashed lines denote the contractions between 
two tensors. 


lattice, which represents the partition function Z, as shown 
in Fig. 2(a). The partition function TN consists of tensors 
7 ’.vi„s 2 ,s3„v 4 located at each vertex, which has four indices .v, 
(i e {1,2,3,4}) corresponding to the four geometric bonds. 
Each index has a finite bond dimension (D = 3), i.e., s, e 
{0,1,2}; Sj = 0(1) means the absence (presence) of a loop 
bond, Si - 2 represents the presence of a dimer, on the spe¬ 
cific edge Si. We properly initiate the tensor T, to make 
sure that a specific lattice site is either occupied by a loop 
[Ti, l.o.o = Ti, o,i,o = Ty o.o.i = To.i.i.o = T 0 ,i,o,i = T 0 ,o,u = 
exp (—v/T)], a monomer [7o,o,o.o = exp (-p/T)], or by a dimer 
[72,0,0.0 = To, 2 , o,o = To, 0 , 2,0 = To,o,o ,2 = exp(-H/T)], and the 
rest elements are zero (forbidden). 

To calculate the thermodynamics of the MDL model, one 
has to efficiently (and accurately) contract the partition func¬ 
tion TN for calculating free energy per site /, energy per site 
e and other thermodynamic quantities. In this work, we define 
the system on two kinds of geometries: for the infinitely large 
2D lattice system, we adopt the iTEBD method [21, 22] for 
accurate contractions; for the cylindrical geometry with finite 


(small) circumferences and infinite length, we diagonalize the 
spectrum of the transfer matrix and evaluate properties with its 
dominating eigensystems. iTEBD was initially proposed for 
efficient simulations of the time evolution and the ground state 
property (through imaginary-time evolution) of ID quantum 
systems, and then generalized to calculate the thermodynam¬ 
ics of 2D classical statistical models [22] and also ID quantum 
lattice models [23]. In our practical simulations, we perform 
the contraction of MPS with transfer MPO until the prescribed 
convergence criterion is reached, say, free energy per site con¬ 
verges to 10~ 14 (almost machine precision). The total number 
of iterations is around 10 3 ~ 4 , depending on the temperatures 
and the physical parameters of the model. The retained bond 
dimension of the boundary MPS D c ~ 100, the convergence 
with D c is always checked, the truncation error is less than 
10 6 at the critical point, and reaches the machine precision 
away from the critical points. 

Monomer-Loop model. — In the partition function Eq. (1), 
if we forbid the dimer occupation (i.e., u —> oo), the model 
is reduced to a monomer-loop model, which can be related to 
the well-known Ising model. For instance, the triangular lat¬ 
tice Ising model can be mapped to a monomer-loop model on 
its dual honeycomb lattice, where the loops are the magnetic 
domain walls separating spins which have opposite orienta¬ 
tions, and the monomers are the topological excitations on 
top of that [15], In our present model, we treat directly the 
monomer-loop picture, and thus can tune monomer energy p 
continuously (with fixed v = 1), making the simulations be¬ 
yond the exact Ising model mapping (also notice that we are 
simulating the loop models on a square lattice, and do not 
allow loop crossovers, which thus in the beginning lacks an 
exact Ising mapping). 

Firstly, we investigate the monomer-loop model with p = 
0, +0.2. The specific heat Cy curve is shown in Fig. 3(a), 
which is computed by taking first-order derivative (versus 
T) of the energy per site. The latter is calculated via e = 
Z*/Z, where Z* is obtained by contracting the TN with one 
T tensor replaced with an impurity tensor T 1 = pTo,o,o,o + 
v (Ti,i,o,o + Ti,o,i,o + Ti,o,o,i + To, 1 , 1,0 + To.i.o.i + To,o,i,i). In 
Fig. 3(a), divergent peaks appear at T c ~ 0.927,1.157,1.39, 
for p - 0.2,0, -0.2, respectively, suggesting the presence of 
second-order phase transitions. 

In addition, to confirm the existence of the phase transi¬ 
tion, we calculate the correlation length £ of the monomer- 
loop model via the following formula 

f=l/ln(^), (2) 

42 

where A\ ( A 2 ) is the largest (second-largest) eigenvalue of the 
transfer matrix M (in case the largest eigenvalues are n-fold 
degenerate, /E is the n +1 largest one). In Fig. 3(b), we observe 
that the correlation length also shows a divergent peak at T c , 
confirming the criticality at the transition point. 

However, interestingly, we find no local order parameters 
for detecting this phase transition, since both the high- and 
low-7’ phases are disordered and have no symmetry breaking. 




















































3 




FIG. 3: (Color online) (a) The specific heat Cy of the monomer- 
loop model with p = 0, ±0.2 and v = 1, where T c ’s correspond to 
0.927,1.157, 1.39 respectively, (b) The correlation length f versus 
temperature T of the classical loop model. The heights of the peaks 
at T c grow by increasing D c . 


We show the numerical results of bond density (n A and hr) 
in Fig. 4 (a), which counts the average bonds (of the loops) 
per site. The results (with p = 0) are shown in Fig. 4(a), from 
which we can see that the low-7’ phase has relatively low bond 
density and thus can be regarded as monomer-condensation 
(MC), while the high-7’ region is a loop-condensation (LC) 
phase. Although n A and hb change from zero to nonzero val¬ 
ues when T increases, they change smoothly through the tran¬ 
sition point. In addition, the same bond densities n A — n B are 
observed for all temperatures T, which suggests that the sym¬ 
metry between two sub-lattices is also intact. Therefore, we 
conclude that the bond density n can not serve as a local order 
parameter for distinguishing two phases. Besides, in Fig. 4(a), 
we also show the bond density n — 1.602 944603 316996 
(with 16 converged significance digits) in the T — °o limit, 
where the state is an equal-weight (classical) superposition of 
all possible monomer-loop configurations. Compared to the 
dimer density n,j = 0.638 123 109 228 547 in monomer-dimer 
model [8], we find that n > 2nd here. 

We also investigate the entropy of the system, including two 
kinds of entropies, i.e., the conventional thermodynamic en¬ 
tropy S — (U - F)/T and the formal “entanglement entropy” 
S e evaluated from the boundary MPS. The latter can be ob¬ 
tained by Se = — E/A^ln(Ar), where A,’s are the Schimidt 
spectrum of the decomposition on a bond i. As shown in 
Fig. 4(b), the bipartite entanglement entropy Se of the clas¬ 
sical loop model shows a clearly divergent peak at T c , indi¬ 
cating the occurrence of a phase transition. On the contrary, 
the conventional thermodynamic entropy S is smooth around 
T c . However, its first-order derivative has a divergent peak, as 
shown in the inset of Fig. 4(b), which is not surprising since 
= -y. Therefore, we see that this “entanglement entropy” 
Se is more sensitive to the phase transition (than the thermo¬ 
dynamic entropy S ), and serves as a very useful numerical 
tool detecting phase transitions. Similar behaviors have al¬ 
ready been seen in our previous tensor-network study of the 
monomer-dimer model [8]. 

In order to further investigate the phase transitions, we also 
define the MDL model on cylindrical geometries. On the 
cylinders with finite circumferences (and infinite length), we 


FIG. 4: (Color online) (a) The mean bond occupation number n A and 
n B on the vertices A and B in the monomer-loop model (p = 0, v = 
1), the horizontal dashed line shows the bond density in the T = oo 
limit. The inset amplifies and shows the low temperature behaviors, 
(b) The entropy S and the entanglement entropy 5 £ of the classical 
loop model with p = 0, v = 1. Inset shows the first-order derivative 
of entropy dS / dT. 

can contract the TN exactly and thus obtain the thermody¬ 
namic quantities. Specifically, we start from both ends, and 
contract the boundary vectors with the transfer matrix M con¬ 
sisting of a column of rank-4 tensors [see Fig. 2(b)]. Repeat¬ 
ing the contraction process until both boundary vectors con¬ 
verge, with which we can evaluate the observables like the en¬ 
ergy expectation values. As the cylinder widths increase, the 
observables should eventually converge to the thermodynamic 
limit results obtained with iTEBD contractions above. 

We calculate the (normalized) gap of the transfer matrix 
6 = | 1 for various cylinder widths w, which is shown in 

Fig. 5(a). In the T > T c region, we observe two degenerate 
dominating eigenstates in the spectrum of transfer-matrix M. 
In particular, as shown in the inset of Fig. 5(a), 6 extrapolates 
to zero at critical point T - T c , in the thermodynamic limit. 
Therefore, 6 can be taken as an order parameter detecting the 
phase transition between the low- and high-T phases. 

We are also interested in the parity of the dominating eigen¬ 
vector x °f transfer matrix M. Since the loops are closed in 
the MDL model, M conserves the parity symmetry. When the 
cylinder is cut into two halves vertically, the number of inter¬ 
sected bonds by the cut is either even or odd, which defines the 
parity of eigenvector^-. For cylinders with open ends (i.e., no 
dangling bonds on the edges), all the allowed configurations 
constitute the even sector, and the dominating eigenvector \e 
in this sector is with even parity. On the other hand, if we 
introduce odd number of open strings on the cylinder, stretch¬ 
ing from the very left boundary to the rightmost side, then 
all the allowed configurations constitute an odd sector, with 
dominating eigenvector;^ of odd parity. Fig. 5(a) shows that 
the dominating even and odd eigenvectors (x e and^v,) become 
degenerate when T > T c . 

Furthermore, we introduce a string operator 0 winding 
around the cylinder to measure the parity of eigenvectors x- 
The string operator 0 = ]~[/Li Pi is a product of the operator 


which lives on the horizontal edge [20]. P — 1(—1) if the 
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FIG. 5: (Color online) Topological characterization of the monomer- 
loop model for various cylinder widths w. (a) The gap of the transfer 
matrix 6 vanishes when T > T c , the inset shows that the extrapolated 
<5 — 0 in the w = oo limit, (b) The odd string order parameter 0 O is 
zero when T <T C and nonzero T > T c ; 0 e also changes its behavior 
at T c . 



FIG. 6: (Color online) The topological properties of the monomer- 
loop model with p = ±0.2 and v = 1 for various cylinder widths 
w. (a) The gap of transfer matrix S vanishes when T > T c , the inset 
shows 6 extrapolates to zero at T c . (b) The string order parameters 
©o and ©,,. 

edge is not occupied (occupied by a bond). Therefore, the 
expectation value of the product of P tells whether the sys¬ 
tem is in the even or odd sector. In Fig. 5(b), we thread an 
open string in the cylinder, and show the numerical results 
Of &o = \Xo®Mx* 0 KXoMxl + XeMx*e) I for various cylinder 
widths, where the partition function Xe(Xo) is the even(odd) 
dominating eigenvectors. We observe that 0 O is a constant 
zero for T < T c , and becomes nonzero when T > T c . In the 
meantime, we also show the even 0 f = [ \ e ®Mx* e / (x<Mx* 0 + 
XeMx* e )\, which is a constant one in the trivial MC phase, 
while also changes its behavior at T c . Thus, © 0 can also be 
taken as an order parameter for detecting the phase transition, 
called nonlocal string order parameter. 

In addition to p = 0 case. Fig. 6 shows the corresponding 
results of the monomer-loop model with p = +0.2, including 
the (normalized) gap 6 in Fig. 6(a) and SOP in Fig. 6(b). Sim¬ 
ilar behaviors can be seen as in the p = 0 case. In Fig. 7, 
we tune the monomer doping parameters p, collect the phase 
transition points of the monomer-loop model with various pa¬ 
rameters p, and obtain the p — T phase diagram. When p < 1, 
there exist second-order phase transitions separating the low- 
and high-7’ phases. T c decreases with increasing monomer 
energy until p = 1, where T c = 0, i.e., no phase transitions. 

Here we would like to address some remarks on the classi¬ 
cal topological order in the monomer-loop model. As a con¬ 



FIG. 7: (Color online) The phase diagram p - T of the monomer- 
loop model. The disordered monomer-condensation (MC) and topo¬ 
logically ordered loop-condensation (LC) phases are separated by a 
second-order transition line when p € (-oo, 1], 

sequence of the loop condensation [Fig. 1(b)], there exist two 
degenerate eigenvectors^ in the high-7’ phase, meaning that 
the phase space is decomposed into two distinct topological 
sectors, which leads to a topological ergodicity breaking. One 
gets exactly the same results by evaluating the thermodynamic 
quantities in either sector, while it is not possible to shift from 
one sector to the other by changing the loop configurations 
only locally. This glassy behavior of the LC phase is due to 
topological reasons, therefore the LC phase can also be called 
a topological glass [15], and the phase transition between MC 
and LC phases is thus a topology-driven transition. 

Monomer-Dimer-Loop model .— In the monomer-loop 
model, the dimer occupation was not allowed (i.e., effectively 
u — oo in Eq. 1). The dimer can be regarded as the “minimal” 
loop of length two [shown in Fig. 1 (c)]. In the following, 
we switch on dimer coverings, and study the full MDL model 
with u — 3 (and p = 0, v = 1). 

As in the monomer-loop model case, we also calculate the 
specific heat Cy, the correlation length <f, bond density n and 
the entropy (the thermodynamic entropy S and the entangle¬ 
ment entropy S e) of the system [20]. CV, S h and the first- 
order derivative of the thermodynamic entropy all show a di¬ 
vergent peak at T c ~ 1.179, indicating the occurrence of a 
second-order phase transition. On the other hand, the bond 
density (the presence of a dimer is considered to be equivalent 
to two overlapping bonds) is smooth around T c and n,\ = n H . 
Thus the bond density again can not serve as a local order 
parameter. 

The normalized gap 6 and SOP ©„ in the MDL model are 
shown in Fig. 8. In Fig. 8(a), we again see a twofold degen¬ 
eracy in the transfer matrix spectrum in the LC phase, indi¬ 
cating that the LC phase is also topologically ordered in the 
MDL model. The inset of Fig. 8(a) shows 6 extrapolates to 
zero (in the thermodynamic limit) at T c . In Fig. 8(b), 0„ ( . are 
shown, in which ©„ is nonvanishing in the LC phase, suggest¬ 
ing that it can also identify the classical topological order in 
the MDL model. In summary, similar to the monomer-loop 
case, the classical topological order exists in the LC phase of 
the general MDL model and the second-order phase transition 
separates the trivial disordered MC and topologically ordered 
LC phases. 6 and 0„ can be used as order parameters to char¬ 
acterize the topology-driven classical phase transition well. 
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FIG. 8: (Color online) The topological properties of the MDL model 
with fi = 0, v = 1, u = 3 for various cylinder widths w. (a) The gap 
of transfer matrix 6 serves as an order parameter; Inset: 6 at T c for 
various cylinder widths and their extrapolation, (b) The string order 
parameter 0 O and ©,,. 

Conclusion and outlook .— In this work, we have system¬ 
atically studied the classical loop model with monomer and 
dimer doping. Using the boundary MPS contraction method, 
we evaluate the partition function TN and obtain the thermo¬ 
dynamic properties including the specific heat Cy, correlation 
length £ and entropies. There exist second-order phase tran¬ 
sitions separating the trivial monomer-condensation and the 
loop-condensation phases, which can not be described by the 
local order parameters like bond density n. However, in the 
LC phase, we find twofold degenerate dominating eigenval¬ 
ues in the transfer matrix spectrum, one in even and the other 
in odd topological sectors. The existence of two topological 
sectors actually breaks the ergodicity. The non-vanishing non¬ 
local order parameter SOP ©„ can also be used to distinguish 
two sectors and thus detect the phase transition. Therefore, 
these two phases can be identified by their distinct topological 
properties, and the phase transition between them belongs to 
a topology-driven type. 

Besides the closed loop cases studied in the MDL model 
above, it is also interesting to consider the model with the 
branching loops [see Fig. 1(d)], i.e., the classical string-net 
model. It is quite straightforward to generalize the tensor- 
network representation here to the classical string-nets, and 
our preliminary calculations show that there also exists a 
second-order phase transition between LC and MC phases. 
However, owing to the existence of branching loops, the 
transfer-matrix breaks the parity symmetry and no longer has 
the well-defined even and odd topological sectors as the MDL 
model has here. Our study of the classical string-nets will be 
published elsewhere. 
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Supplementary materials for “Topology-driven phase transitions in the classical 

monomer-dimer-loop model” 


TENSOR NETWORK REPRESENTATION 


In this part, we introduce the partition function TN representation of these models and the contraction method of the nonlocal 
string order parameter 0 D , e . 

Since each vertex is allowed to be covered by at most one loop or a “thick” dimer in the MDL model, there are eleven 
nonzero elements in the vertex tensor T. The allowed nonzero tensor elements T Su s 2 , i3 , S4 and their corresponding classical 
loop configurations are schematically shown in Fig. SI (a-k). Tq,o,o,o — exp (-/3/j) corresponds to the absence of any loop, 
i.e., a monomer [Fig. SI (a)], T i,i,o,o = 2”t.o.i.o = 2”t.o.o.i = 7o,i,i,o = 7o,i,o,t = 7o,o,i,i = exp(-jSv) represent the one-loop 
configurations [Fig. SI (b-g)], 72 ,i,i,i = 7”i, 2 ,i,i = T 14 , 2,1 = T 1 , 1 , 1,2 = exp(-/?«) describe the vertex with a dimer [Fig. SI (h-k)], 
and Ti 110 = 7o,i,i,i = 74, 0 , 1,1 = 74, 1 , 0,1 = exp(-ySv) introduce the branching loop configurations [Fig. SI (l-o)]. 

In Fig. S2, we show the way of evaluating 0„ f; on a cylindrical geometry. Denominator is the partition function, and the 
numerator measures the expectation value of 0 in the odd(even) sector. 


SPECIFIC HEAT, CORRELATION LENGTH, AND ENTROPY RESULTS OF A MONOMER-DIMER-LOOP MODEL 

For the MDL model with /r = 0, v = 1, u = 3, we also calculate its thermodynamic properties: the specific heat Cy, the corre¬ 
lation length £, the bond density n and the entropy that includes the traditional thermodynamic entropy S and the entanglement 
entropy Se- In this part, we show the results as follows. 

Fig. S3 shows our computed specific heat Cy of the MDL model with the energy u — 3 of a dimer and fixed \x = 0, v = 1, 
where a divergent peak of Cy occurs at T c ~ 1.179, uncovering the existence of the second-order phase transition. 

In Fig. S4, we show the correlation length of the MDL model. From Fig. S4, we can also observe a divergent peak at T c , the 
second-order phase transition point. 

The bond densities per site ua and n b on vertexes A and B in the MDL model with p = 0, v = 1, h = 3 are shown in Fig. S5. 
Both ha and 11 h are changing smoothly around T, and n A is equal to 11 /-; for all temperatures, meaning that the symmetry between 
two sub-lattices in the system is not broken. 

The thermodynamic entropy S and the entanglement entropy S e of the MDL model with /x — 0, v = 1, u = 3 are shown 
in Fig. S 6 . S e shows a divergent peak at T c , while S is smooth around T c , and its singularity can only be seen after taking a 
first-order derivative over T, which is shown in the inset of Fig. S 6 . 
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FIG. SI: (Color online) The eleven (allowed) nonzero elements of the vertex tensor T SUS2jS}tS4 in the partition function TN representation and 
their corresponding loop configurations. In the TN, a lattice site occupied by: (a) a monomer; (b-g) a loop; (h-k) a “thick” dimer; (l-o) a 
branching loop. 
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FIG. S2: (Color online) The contraction method of the nonlocal string order parameter 0 oe on a cylindrical geometry, corresponding to the 
equations for 0 D e in the main text. The definition of x„,e> M, 0 is also explained in the main text. 
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FIG. S3: (Color online) The specific heat Cy of the MDL model with the energy u = 3 of a dimer and fixed fi = 0, v = 1. 
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FIG. S4: (Color online) The correlation length of the MDL model with /j = 0, v = \,u = 3. The heights of the peaks at T, grow with the 
increase of D c . 



FIG. S5: (Color online) The bond density per site n A and n B of the MDL model with ^ = 0, v = 1, u = 3. Inset: The amplification of n A and n B 
at low temperature region. 



FIG. S6: (Color online) The thermodynamic entropy S and the entanglement entropy Se of the MDL model with /t = 0, v = l,u = 3. The 
inset shows the first-order derivative of entropy S for temperature dS/dT. 





















